The Elevation and Duration of Wave Crests 


WHO] 


by 
William N. Seelig, John P. Ahreng§, DOCUMENT 
and COLLECTION 


William G. Grosskopf 


MISCELLANEOUS REPORT NO. 83-1 


JANUARY 1983 


U.S. ARMY, CORPS OF ENGINEERS 
COASTAL ENGINEERING 
RESEARCH CENTER 


Kingman Building 
ach. Fort Belvoir, Va. 22060 


Reprint or republication of any of this material 
shall give appropriate credit to the U.S. Army Coastal 


Engineering Research Center. 


Limited free distribution within the United States 
of single copies of this publication has been made by 
this Center. Additional copies are available from: 

Nattonal Technical Information Service 
ATTN: Operattons Diviston 
5285 Port Royal Road 
Springfield, Virginia 22161 
Contents of this report are not to be used for 
publication, or promotional purposes. 


advertising, 
Citation of trade names does not constitute an official 
of the use of such commercial 


endorsement or approval 
products. 


The findings in this report are not to be construed 
as an official Department of the Army position unless 


so designated by other authorized documents. 
: | 


1 


| 


NOILIIT1I09 
INAWND0G 


l 


COT 


ni 


uN 


UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


READ INSTRUCTIONS 
REPORT DOCUMENTATION PAGE BEFORE COMPLETING FORM 
T. REPORT NUMBER 2. GOVT ACCESSION NO, 3. RECIPIENT'S CATALOG NUMBER 
MR 83-1 


4. TITLE (and Subtitle) 5. TYPE OF REPORT & PERIOD COVERED 


Miscellaneous Report 


THE ELEVATION AND DURATION OF WAVE CRESTS 


6. PERFORMING ORG. REPORT NUMBER 


8. CONTRACT OR GRANT NUMBER(s) 


7. AUTHOR(s) 
William N. Seelig 
John P. Ahrens 
William G. Grosskopf 
9. PERFORMING ORGANIZATION NAME AND ADDRESS 
Department of the Army 

Coastal Engineering Research Center (CERRE-CS) 
Kingman Building, Fort Belvoir, VA 22060 
11. CONTROLLING OFFICE NAME AND ADDRESS 


10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 


A31592 


12. REPORT DATE 
January 1983 
13. NUMBER OF PAGES 


73 


1S. SECURITY CLASS. (of thie report) 


Department of the Army 
Coastal Engineering Research Center 


14. 


UNCLASSIFIED 
1Sa. DECL ASSIFICATION/ DOWNGRADING 
SCHEDULE 


DISTRIBUTION STATEMENT (of thie Report) 


Approved for public release; distribution unlimited. 


DISTRIBUTION STATEMENT (of the abstract entered in Block 20, if different from Report) 


» SUPPLEMENTARY NOTES 


- KEY WORDS (Continue on reverse side if necessary and identify by block number) 


Wave breaking 
Wave crests 


Irregular waves 
Monochromatic waves 


ABSTRACT (Continue em reverse side if neceasary and identify by block number) 
The stream-function wave theory of Dean (1974) is used together with 
monochromatic and irregular laboratory wave data to develop methods for 
estimating the elevation and duration of wave crests. The resulting predic- 
tion techniques are applied to a wide range of wave conditions measured at 
CERC's Field Research Facility in Duck, North Carolina, and are shown to give 
reliable and often conservative estimates of crest elevation. The techniques 
presented in this report can be used for both nonbreaking and breaking wave 
conditions. 


FORM 
DD . jan 7a 1473 EDITION OF 1 Nov 65 1s OBSOLETE 


UNCLASSIFIED 
SECURITY CLASSIFICATION GF THIS PAGE (When Data Entered) 


PREFACE 


As waves move toward the shoreline they become increasingly nonlinear with 
the elevation of the wave crests becoming relatively large and short in dura- 
tion. This report is published to provide coastal engineers with techniques 
for estimating crest elevation and duration. The techniques were developed 
using laboratory data from a number of sources covering a wide range of con- 
ditions for both monochromatic and irregular waves under nonbreaking and 
breaking conditions. Predictions compare favorably with wave observations 
made at the U.S. Army Coastal Engineering Research Center's (CERC) Field 
Research Facility at Duck, North Carolina. The work was carried out under 
CERC's Wave Estimation for Design work unit, Coastal Flooding and Storm Pro- 
tection Program, Coastal Engineering Area of Civil Works Research and Develop- 


ment. 


This report was prepared by William N. Seelig, Research Hydraulic Engineer, 
John P. Ahrens, Oceanographer, and William G. Grosskopf, Hydraulic Engineer, 
under the general supervision of Dr. R.M. Sorenson, Chief, Coastal Processes 
and Structures Branch and Mr. R.P. Savage, Chief, Research Division. 


Technical Director of CERC was Dr. Robert W. Whalin, P.E. 


Approved for publication in accordance with Public Law 166, 79th Congress, 
approved 31 July 1945, as supplemented by Public Law 172, 88th Congress, 
approved 7 November 1963. 


TED E. BISHOP 


Colonel, Corps of Engineers 
Commander and Director 


CONTENTS 


CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC. . 


SH AVOSKOVESS ZNNI0) IDDVMOMCIILONS 5 6 6 6 6 0.0 00000600 06 
i IORAG NOMI GIIOM 6 6 G6 Oo 0 60 0 6.0 00 og oo Oo ne 
Il MAICHBVAIEONIS TONOAY 9 Gp Gg 19 0, 6 G6 OO 0 OG 0 0 60.6 G6 Oo 9 
IIL SOW OF IWYYNS Go 6) olo 06 6 6 06 6 06:06 00 050 0 
1. Laboratory Tests Conducted by Authors ....... 
2 Other, Kaboratony mlests Usediey. cucu siren tele oie 0 

35 PROroOmoe Waka Use 6 6616 06605560000 


IV TEST RESULTS AND PREDICTION TECHNIQUES FOR MONOCHROMATIC 


1. Monochromatic Wave Crest Elevations... . 
2. Crest Elevation Prediction Aids ..... 


3. Relative Crest Duration for Monochromatic Waves . . 


V TEST RESULTS AND PREDICTION TECHNIQUES FOR IRREGULAR 
1. Irregular Wave Crest Elevations ........ 
2. Crest Elevation Prediction Aids ........ 
3. The Duration of Irregular Wave Crests... . 
VI COMPARISON OF RESULTS FOR PROTOTYPE DATA... . . 
io isteilel Comebticitemss oc 5 6 6 6600000 0 0 0-0 
he LIGNE ANS, WDiesleskeetenteye 5 6 6 56 0 6 6 0 G6 OO 
VIL SUMMARY SAND CONGIEUSTIONS Ie siete circled Welinenn sts to KeMlis 6 
Ug LAV AVEO (GMIVAD Gg o 6. 0 O10 .0'6 6 050.040 6 0 
APPENDIX 
A MONOCHROMATIC LESTSARESUIERS acm vl eoitcniennciircnn stir. : 
B IRREGULAR WAVE CREST ELEVATIONS ...... .- 
C IRREGULAR WAVE CREST DURATIONS ........ © 
D IRREGULAR WAVE DATA FOR SHOALING CONDITIONS . . O 
iE; (EO ICOVNE, VPI, ANOS ISIS Gg og 6 G6 6b 6 6 6 0.0 . 
F ESTIMATING THE WAVE CREST ELEVATION AT BREAKING 
ON: Av PANE. SHORE 9/0) ces glean ui ay ct rates came eo ear rsiaaes 
TABLES 
1 Summary of CERC laboratory test data .....-+-+-e+.-s 


2 Values of the empirical coefficients Ay and A3 for 
GULBIS CIEABIES o 6150 0 6060560060000 00 ¢ 


WAVES .. 


e e e 


predicting 


Page 


50 


Dil 
56 
60 
64 


72 


73 


15 


20 


N DD WN ££ 


10 


11 


fF WH NH Ff 


oOo ON DN 


10 
11 
12 
13 
14 
15 
16 
17 
18 


CONTENTS 
TABLES-—-Continued 


Observed versus predicted relative crest elevations for 


mYepVOOlaaomayesle WEES 6 6°66 5016 06 6 6 6.06006 4 00 0.6 
Coefficients for estimating the duration of wave crests... . 
Observed versus predicted relative crest durations ..... .- 
Selected irregular wave crest elevation data ...... . -» 


Selected irregular wave crest data for waves breaking on a 
i om SO Bil@pS 66 016 40°00 6059 6 6 6 6 6°05 50 0 Oo 


Values of A* for levels of probability of exceedance for 
irregular waves derived from 122 laboratory experiments .. . 


Values of A* for levels of probability of exceedance for 
irregular waves on a 1 on 30 laboratory slope. . ...... «» 


Summary of prototype wave conditions. . .... +. .++ ++. -» 


Values of crest heights observed and predicted for waves 


ale Wels, Nore Carrollsineig:o 56 696 6 Oo 6 6 0650005005 56 0.0 
FIGURES 
Example of an extremely nonlinear wave condition .... . 
Predicted monochromatic wave profile. . ... +. +. +6 «© »« » «= « « 
Experimental setup in CERC's 4.5-meter-wide tank ... . 


Sample of water level time histories taken at four gages during 


Aim sliceGetllehe TENE TESS 6 56 6 46.6 06 6 6060605666 
Experimental setup in CERC's 0.9-meter-wide tank ...... 
Experimental setup in CERC's 0.4-meter-wide tank ...... 
Observed and predicted wave spectra . ... +... e+ ++ « « 


Empirical coefficients A, and A3 for estimating n,/H 
Effect of beach slope on relative wave crest elevation .... 
Effect of beach slope on breaker crest elevation. ..... .» 


Wave crest elevation for constant values of wave steepness . 


Wave crest elevation for constant values of H/Hy, ¢ 0 0 0 
n/d for constant values of wave height to depth ..... 
Empirical parameters for estimating wave crest durations .. . 
Values of dimensionless wave crest duration ........ - 
Sample of irregular wave time histories . ......++. +e « 


Variation in wave crest elevation for al on 30 slope... 


Sample of water level time histories for irregular waves at 
selected locations along a 1 on 30 slope .........-. - 


5 


Page 


20 
28 
28 
31 


33 
38 


38 
46 


46 


10 
11 
ale 


13 
14 
5) 
16 
19 
Zadl, 
22 
23 
25 
26 
27 
29 
30 
33) 


34 


19 
20 
21 
22 


23 
24 


25 


26 


27 
28 
29 


CONTENTS 


FIGURES--Cont inued 


Sample of irregular wave crest elevation distributions 


Sample of wave crest distributions. ..... . 


Values of AT for probability of exceedance for irregular waves 


Observed and predicted crest elevations exceeded by 2 percent 


Of stheliwaviewmCrestsivemyc lenin veuier ienmiern-cii vena chmtelnnonis 


Wave crest elevation exceeded by 13.5 percent of the wave crests 


Crest elevation at the 2-percent probability level of exceedance 


Ope sleeve VERS G25 oo Gio. 6 0 oligo 6 


Crest elevation at the 13.5-percent probability 
exceedance for irregular waves ........ 


Wave crest elevations along a 1 on 30 slope for 
percent probability of exceedance levels... 


level of 


° e 


2- and 13.5- 


Dimensionless wave crest duration for irregular waves . 


Observed versus predicted crest elevations (p = 


Observed versus predicted crest elevations (p = 


0.135) 
0.020) 


CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (SI) UNITS OF MEASUREMENT 


Doe customary units of measurement used in this report can be converted to 
metric (SI) units as follows: 
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THE ELEVATION AND DURATION OF WAVE CRESTS 


by 
William N. Seeltg, John P. Ahrens, and Willtam G. Grosskopf 


I. INTRODUCTION 


Coastal structures are generally built to withstand extreme waves that 
occur infrequently. One of the important aspects of these waves is that they 
may be highly nonlinear if the wave steepness is large or the relative water 
depth is small. Two of the important engineering characteristics of nonlinear 
waves are that the crest elevations are more than half of the wave heignt and 
the duration of the crests are less than half of the wave period. 


Crest elevation is important because the forces on a structure may be more 
strongly related to crest elevation above the water level than the trough dis- 
tance below the water level. Uplift forces on piers and fixed structures 
are particularly sensitive to wave crest elevation, which is one of the reasons 
why most offshore platforms and piers are built well above the water level. 
Crest elevation also influences wave runup and overtopping (Weggel, 1976a). 


The relatively short duration of nonlinear wave crests is important, 
because the combination of high elevation and short duration produces large 
water particle motions and accelerations. These hydraulic characteristics have 
the potential for suspending sediment, lifting armor units, and causing surpris-— 
ingly high drag forces on coastal structures. 


Figure 1 is a sample wave record made as part of the Coastal Engineering 
Research Center's (CERC) coastal wave collection program (Thompson, 1977). 
This example clearly shows the high crest elevations and short crest durations 
that may occur in coastal waters. 


This study develops methods of predicting crest elevations and duration for 
monochromatic and irregular waves. Laboratory waves are used to minimize the 
influence of refraction, diffraction, irregular topography, directional spread 
of wave energy, and multiple wave trains on the characteristics of wave crests. 
Assumed independent variables used in the prediction methods are water depth, 
wave period or period of peak energy density for irregular waves, wave height 
(significant wave height for irregular waves), and beach slope. Predictions 
are compared to selected prototype wave measurements taken at the CERC Field 
Research Facility (FRF) at Duck, North Carolina. 


II. LITERATURE REVIEW 


Goda (1964) performed some basic research on the elevation of monochromatic 
laboratory wave crests. Using his own data and the data of Bretschneider (1958), 
Goda found that the ratio of wave height to water depth played a dominant role 
in determining the relative crest elevation. However, there are considerable 
scatter in the results. 


Jahns and Wheeler (1973) and others have performed theoretical derivations 
of expected crest heights for irregular wave conditions. Without including 
nonlinear effects they conclude that a Rayleigh-type distribution should 
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Figure 1. Example of an extremely nonlinear wave condition 
(recorded by a CERC continuous-wire staff at Lake 
Worth, Florida, 28 March 1971 at 0820). 


provide an upper limit or conservative estimate of crest elevations correspond- 
ing to various probabilities of exceedance. Actual irregular wave data have 
crests that are higher than given by this theory, so Jahns and Wheeler (1973) 
suggest an empirical correction factor. One limitation of the correction 
factor is that the observed probability of exceedance may deviate several or 
many orders of magnitude from the Rayleigh distribution. 


Dean (1974) developed one of the most comprehensive theories to date for 
predicting properties of monochromatic waves traveling over a flat bottom. 
This higher order stream-function theory can be used to predict crest elevation 
and duration for a wide range of wave and water level conditions. The theory 
indicates that in the Airy limit, where the wave amplitude approaches zero, 
the wave crest elevation is approximately equal to one-half the wave height 
and the duration of the crest is half the wave period. In the cnoidal wave 
limit the crest elevation approaches the wave height and the duration of the 
crest becomes small compared to the wave period. Figure 2 is an example 
predicted water level time history for a highly nonlinear wave. In this 
extreme condition the ratio of crest elevation, Nc» to wave height, H, is 
Ne/H = 0.90. The duration of the crest, T is one-fourth the wave period, 
ie he Wt OeL5c 


(or 


Singamsetti and Wind (1980) presented wave crest characteristics for 
monochromatic laboratory wave conditions at the point of breaking for constant 
beach silopes, m, of 1 on) 405 1) on) 205, 1) on 10, vand 1lyon)5). 


ne) 


1-/H = 0.90 d/gT? = 0.0008 
Te/T = 0.25 H/gT* =0.0003! 


0 | Tati. T 
Time —= 


Figure 2. Predicted monochromatic wave profile (after Dean, 1974). 


III. SOURCES OF DATA 


Sources of data upon which the predictive techniques are based include 
laboratory data collected by the authors for a wide range of conditions, 
laboratory data collected by other authors, and field data from the FRF at 
Duck, North Carolina. Conditions under which these data were collected and 
methods of analysis are described below. 


1. Laboratory Tests Conducted by Authors. 
a. Wave Conditions, Data, and Test Setups. Most of the monochromatic wave 


data and all the irregular wave data collected over a horizontal bottom were 
obtained in a 105.7-centimeter-wide channel within CERC's 1.2-meter-high by 
4,5-meter-wide by 42.6-meter-long wave tank. The monochromatic wave data are 
tabulated in Appendix A and the irregular wave data are tabulated in Appendixes 
Band C. A gravel absorber beach with a slope 1 on 12 was used at the end of 
the channel to reduce wave reflection to a low level. Four paralleled wire 
resistance-type wave gages were used in the channel to record the wave condi- 
tions. A plan view and a profile view of the test setup with gage locations 
are shown in Figure 3. Figure 4 shows a short section of typical irregular 
wave records taken approximately simultaneously at each of the four gages. 
The time axis of the records has been shifted slightly to place the highest 
crest in the center of each section to show the differences in the time 


It 


Tank Woll 


Test Channel oN 


cole Other Test E 
Wove Goges Channels 


wove Tank Station (m) 19m to 
Wove Generator 
a. Plan View 
Wave sabe Incident Waves 
@ ©® 


____20.4m to Wave 
Generator Blade 


lon l2 Slope 
30cm feral 90 em—ater>} 30 em 


Grovel Absorber 
b. Profile View 


Figure 3. Experimental setup in CERC's 4.5-meter-wide tank. 
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Sample of water level time histories taken at four gages during 
an irregular wave test (time axis was shifted for the highest 
crest to appear at the center of each plot). 


sequences of the records. Since gages 1 and 4 are the farthest apart they 
show the greatest differences. The data from gages 1 and 4 are tabulated in 
Appendixes A, B, and C. 


Some additional monochromatic wave data over a horizontal bottom were 
collected in CERC's 0.9-meter-high by 0.9-meter-wide by 45.7-meter-long wave 
tank. This tank has a 1 on 30 hogs hair absorber beach to keep reflected wave 
energy to very low levels. Three parallel wire resistance-type wave gages were 
used to record the wave conditions. A profile view of the test setup is shown 
in Figure 5; the data are tabulated in Appendix A. 
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Figure 5. Experimental setup in CERC's 0.9-meter-wide tank. 


Wave data were also collected in CERC's 0.9-meter-high by 0.4-meter-wide 
by 45.7-meter-long wave tank. The setup in this tank was designed to obtain 
data on the shoaling and breaking of irregular waves. Ten parallel wire wave 
gages were used: three gages over the horizontal tank bottom and seven over a 
1 on 3 concrete slope. Figure 6 shows a profile view of the test setup and 
gage locations; the data are tabulated in Appendix D. 


Table 1 provides a summary of the wave conditions, data and test setups 
used in the laboratory tests. 


b. Procedures and Analysis. All the waves were generated using hydrauli- 
cally actuated piston-type wave makers. The monochromatic waves studied 
included waves with sinusoidal blade motion with Ursell numbers less than 25 
and cnoidal waves with Ursell numbers greater than 25. The Ursell number, Ups 
is defined 


HL 
Up = Hha (1) 


where H is the average wave height over the flat bottom part of the wave 
tank and Lag the local wavelength, calculated using linear theory and defined 


by 
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Figure 6. Experimental setup in CERC's 0.4-meter-wide tank 


Table 1. Summary of CERC laboratory test data. 


Wave tank Wave 
conditions 


1.2 by 4.5 by 42.6 m | Monochromatic | Horizontal 


1.2 by 4.5 by 42.6 m | Irregular Horizontal 
0.9 by 0.9 by 45.7 m | Monochromatic | Horizontal 


0.9 by 0.4 by 45.7 m Irregular Horizontal 
and 1 on 30 


T? 2nd 
Ly = a tanh i (2) 


where T is the wave period and d the water depth. For sinusoidal blade 
motion the restriction on the Ursell number was to eliminate the influence of 
secondary waves; the cnoidal wave data were collected by Weggel (1976b) and 
are discussed in Section 2. 


The irregular waves used in this study included a wide range of spectral 
types and relative depths. Spectra were simulated by taking various theoreti- 
cal spectra, such as Pierson-Moskowitz, Joint North Sea Wave Project (JONSWAP), 
Bretschneider, and the six-parameter spectra of Ochi and Hubble (1976), and 
slicing the spectra into 60 equal area segments with each segment represented 
by a sinusoidal component at the frequency of the midpoint of the segment. 
Each component was transformed into a piston blade motion using Biesel's 
equation (Biesel, 1951) at the given frequency and assuming a random-phase 
relation among the components. The time histories of the 60 components were 


1S 


linearly superimposed to determine the generator blade motion. Figure 7 shows 
a sample of observed and predicted wave spectra. 
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Figure 7. Observed and predicted wave spectra (sample of 
Pierson-Moskowitz spectra). 


In all tests the water surface time history was recorded at each wave gage 
for 4,096 data points. A sampling rate of 16 times per second was used for 
waves with periods less than 3 seconds, and a rate of 8 times per second was 
used for waves with periods greater than 3 seconds. The data were digitized 
and recorded on magnetic tape through the use of a Data Acquisition System 
(DAS). An analysis of the tapes was performed later on a general-purpose 
computer. 


Two distinct types of data analysis were performed on the digitized wave 
records: a fast Fourier transform (FFT) and a zero up-crossing analysis. The 
FFT was used to determine the spectrum of the wave record, such as shown in 
Figure 7, which in turn was used to determine the period of peak energy 
density. The period of peak energy density, Tp, is the reciprocal of the 
frequency of the midpoint of the 11 consecutive spectral lines with the most 
wave energy. The zero up-crossing method of analysis is used to determine the 
wave heights, crest heights, wave periods, and crest durations of individual 
waves in irregular wave trains or average values for monochromatic wave trains. 
Crest height, n., is measured relative to the mean water level (MWL) and 


16 


crest duration, T,, is the time difference between an adjacent zero up-cross 
and a zero down-cross; zero indicates the MWL of the record. The significant 
wave height, H,, for irregular waves was computed by multiplying the root- 
mean-square (rms) water surface displacement by four. All the analysis was 
based on the assumption that there was no wave reflection in any of the 
channels where data were collected. 


For irregular wave conditions the elevations of the wave crests were 
ranked from the highest to the lowest and the probability of exceedance, p, 
for various levels of probability were calculated. For example, N9.01 indi- 
cates a crest height with a probability of being exceeded by 1 percent of the 
crest heights, i.e., p= 0.01. Since the larger waves were of the most 
interest the following probability levels for crest heights were calculated: 

p = 0.005, 0.01, 0.02, 0.05, 0.10, 0.135, 0.15. Dimensionless wave crest 
height probability levels, (nc)p/Hs, for wave gages 1 and 4 of the test setup 
shown in Figure 3 are tabulated in Appendix B. The mean dimensionless crest 
duration, (T./T),, was also determined for the highest 5, 10, 15, 20, 25, 33, 
and 100 percent oe the wave crests for gages 1 and 4 and is tabulated in 
Appendix C. 


2. Other Laboratory Tests Used. 


Weggel (1976b) conducted tests using nonsinusoidal blade motion to generate 
cnoidal waves. These waves have Ursell numbers between 27 and 184 and are 
relatively free of secondary waves. The cnoidal wave conditions are given in 
Appendix E. 


Singamsetti and Wind (1980) studied the breaking of monochromatic waves 
on a number of smooth laboratory slopes. The relative crest elevations at 
breaking Tp / Hp and the relative crest durations at breaking, GRyit). are 
tabulated in Appendix A, where nm, is the height of the crest above MWL at 
breaking, H, is the wave height at breaking, and T is the period of a 
monochromatic wave. 


3. Prototype Data Used. 


Prototype wave data collected at FRF, Duck, North Carolina, are used in 
this study, particularly the data measured during a major storm in October 1980 
in which a large variety of wave steepnesses were observed. The gages were 
Waverider buoys located 0.5, 3, 6, and 12 kilometers from shore in an area 
where the bottom contours were generally straight and parallel. The dominant 
wave direction was approximately along the line of the instruments. The data 
were analyzed in the same manner as that previously described for the labora- 
tory data. 


IV. TEST RESULTS AND PREDICTION TECHNIQUES FOR MONOCHROMATIC WAVES 


Monochromatic laboratory wave test results are compared with Dean's stream- 
function wave theory (1974) and the resulting prediction techniques are 
discussed. 


1. Monochromatic Wave Crest Elevations. 


a. Waves Traveling Over a Horizontal Bottom. Dean (1974) developed a 


stream-function wave theory to predict a number of wave characteristics for 
waves propagating over a flat horizontal bottom for a number of relative depths 
and wave heights equal to 0.25, 0.50, 0.75, and 1.00 of the breaking limit. 

An examination of the tabular values of n,/H (Dean, 1974, Vol. II) shows that 
relative crest height can be approximated using the empirical relation: 


Ne a Al (3) 


He lel aniad $e) (aki Rina el 
1.0 + tanh |-A> fn (. x) 
fo) 


where Lo is the deepwater wavelength given by Airy theory as 


2 
ye (4) 


1) Dir 


Wave trough elevation, mm, is equal to wave height minus crest elevation. 
Dean's tabular results are compared with equation (3) and A, is found to be 
1.0 for wave traveling over a flat bottom (m = 0.0). Values of A» and A3 
are functions of relative depth (d/L, or d/(gT?)), as shown in Figure 8 and 
given in Table 2. 


Laboratory data in this study show excellent agreement with Dean's stream- 
function theory and give a value A, = 0.992 + 0.0393 (Table 3). Cnoidal waves 
(Up > 25) have slightly smaller values of this empirical parameter with Aj = 
0.95 + 0.038 (Table 3). 


b. Waves on a Slope. Singamsetti and Wind's (1980) data are used to 
evaluate A, for monochromatic waves at the breaking point on various beach 
slopes. Wave crest elevations, as indicated by A,, decrease approximately 
linearly with beach slope, m (Table 3, Fig. 9). For example, waves at the 
breaking point on a 1 on 5 slope (m = 0.20) have only 67 percent of the crest 
elevations of a wave with same H/ (gT2) and d/(gT) traveling over a flat 
bottom. 


Figure 10 illustrates the predicted influence of beach slope on the 
breaking crest elevation for a sample wave tank condition. In this case waves 
in the incident flat part of the tank have d/(gT?) = 0.019 and H/d = 0.37. 
The figure shows the predicted crest elevation at the breaking point using the 
method described in Appendix F for estimating Ng: In this example slopes 
flatter than 1 on 20 have little influence on the crest elevation, while slopes 
steeper than 1 on 20 produce lower crest elevations. 


2. Crest Elevation Prediction Aids. 


Equation (3), which is used for predicting monochromatic wave crest eleva- 
tions, can be presented in a number of graphical forms useful for predicting 
and understanding trends in relative crest elevation. Figure 11 shows Ne /H 
versus relative depth for various wave steepnesses. Relative crest elevation, 
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Table 2. Values of the empirical coefficients 
Ag and A3 for predicting relative 
crest. 1} 


0.000796 
0.00159 
0.00318 


0.00796 
0.0159 
0.0318 
0.0796 
0.159 


lElevation for horizontal bottoms, using 
equation (3) and A; = 1.0. 


Table 3. Observed versus predicted relative crest elevations for 
monochromatic waves. 


Type of Ay Std. Coefficient Up Source 
Wave dev. of range 

generator variation ae 

motion (pet) 


Sine (9.992 | 0.0393 This study 


Cnoidal 00 0.95 0.038 4.0 13 Weggel (1976b) 


Sine 40 | 0.909! | 0.0375 (Aga 26 Singamsetti 
and Wind 
(1980) 


Sine 20 | 0.88! | 0.035 4.0 Ww Singamsetti 
and Wind 
(1980) 

Sine 10 | 0.818! | 0.0752 9.2 22 Singamsetti 
and Wind 
(1980) 

Sine 5 |) Oo672~ |) OsOs7/0) 365 21 Singamsetti 
and Wind 
(1980) 


lat the breaking point. 
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Figure 9. Effect of beach slope on relative wave crest elevation. 
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Wave crest elevation for constant values of wave steepness. 


Figure 11. 


n./H, is strongly influenced by d/(gT*) for a given value of wave steepness 
for d/(gT*) > 0.05. The breaking limit, H,, given in this figure is from 
Goda (1975): 


d pee 
Hy, = 0.17 Lo {1.0 - exp (-4.712 = [1.0 + 15m (5) 
Oo 


Figure 12 presents n,/H for constant values of the ratio of wave height 
to breaker height, H/H,. The figure shows that in relatively deep water, 
No /H only deviates from the Airy condition of Ne /H = 0.5 when the wave 
height becomes a significant fraction of the breaking wave height. Relatively 
shallower water causes a wave to become more nonlinear, even when the wave 
height is small compared to the maximum possible breaker height. For example, 
ne/H = 0275 Vat d/i(et) =/0,0015s stor) H/o =iOdy(raice'2)e 


The parameter n,/d is presented in Figure 13 for constant values of H/d 
at various values of d/(gT*). This figure shows the combined influence of 
wave steepness and water depth on crest elevation. n./d reaches a minimum 
for constant values of H/d at d/(gT*) in the neighborhood of 0.03. At 
values of d/(gT) larger than 0.03 the increasing wave steepness causes n./d 
to increase until the breaking point is reached. For smaller values of 
d/(gT?) the decreased depth or increased period produces greater wave non- 
linearity and increased n,/d. 


3. Relative Crest Duration for Monochromatic Waves. 


As waves become increasingly nonlinear the duration of the crest, T¢,, 
decreases to less than half of the wave period. Values of relative crest 
duration, Wallaby predicted by stream-function theory (Dean, 1974) can be 
approximated by the empirical relation: 


qu | Lo 
ae = Cj, tanh E Ln (c; =)| (6) 


where C,, Cj, and C3 are empirical coefficients with C, = 0.5 for waves on a 
flat bottom (m = 0.0). Values of Cy and C3 are given in Figure 14 and 
Table 4. Good correlation is found between equation (6) and experimental data 
for small beach slopes as shown in Table 5. However, the equation does not 
apply to beach slopes steeper than 1 on 20 (m = 0.05). 


Figure 15 presents predicted values of T,/T as functions of d/(gT?) 
and H/(gT*). Note that the duration of the wave trough, T; is given by 


(7) 


V. TEST RESULTS AND PREDICTION TECHNIQUES FOR IRREGULAR WAVES 


Irregular wave crest elevations and durations are more variable than mono- 
chromatic wave conditions because wave energy is distributed over a range of 
frequencies. This energy distribution produces waves with a variety of heights 
and periods, which have varying amounts of nonlinearity. Larger waves are of 
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Figure 13. Ne/d for constant values of wave height to depth. 
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Empirical parameters for estimating wave crest durations (m = 0.0). 


Figure 14. 


Table 4. Coefficients for estimating the 
duration of wave crests. 


QO. 9.00032 | 
0.00080 
0.0016 
0.0032 


0.0080 
0.016 
0.032 
0.080 
0.16 
0.32 


relent ere Sie) San IS)h SSS NS) 


1c, = 0.5 for waves on a flat bottom 
(m = 0.0). 


Table 5. Observed versus predicted relative crest 
durations (monochromatic waves-sinusoidal 
generator blade motion). 


Coefficient JO. Source 
of 
variation 
of Cy 
(pet) 


This study 


Singamsetti 
and Wind 
(1980) 


Singamsetti 
and Wind 
(1980) 


Singamsetti 
and Wind 
(1980) 


Singamsetti 
and Wind 
(1980) 


lat the breaking point. 
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Values of dimensionless wave crest duration. 


Figure 15. 


greatest interest, so the primary emphasis of this analysis is on the highest 
15 percent of wave crests in a wave train. The variable p is used to desig- 


nate the probability level of exceedance of interest. For example, 1.) p=0 05? 


where the subscript p = 0.05, indicates that 5 percent of the crest elevations 
exceed the value of Me, 


Assumed independent variables are significant wave height, H, (defined as 
four times the rms water level, n,,,), period of peak energy density, Tp» 
water depth, d, and beach slope, nm. 


di. Irregular Wave Crest Elevations. 


Figure 16 presents several irregular wave records illustrating various 
amounts of nonlinearity. Records A and B have crests which are slightly larger 
than half of their zero up-crossing wave heights; record C has crest elevations 
which are a large fraction of wave height. Table 6, which is a compilation of 
selected wave record results, indicates the elevation of maximum crests for 
various water depth, significant wave height, and period of peak energy density 


d/gTp* = 0.025 Hs/gTp* = 0.0063 H,/d=0.25 
A 

d/gTp* = 0.017 Hs/gTp* = 0.0035 Hg /d = 0.21 
B 

d/gTp* = 0.004 Hs/gTp* = 0.00106 H,/d= 0.26 


Figure 16. Sample of irregular wave time histories. 
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combinations. Thompson (1980) discussed several coastal gage records with large 
ratios of H,/d that have n,/H, ratios on the order of one for the highest 
crest observed, while No/d for the largest crest is somewhat less than one. 
Jahns and Wheeler (1973) found similar extremes for hurricane and storm waves 

in the Gulf of Mexico. Large values of No /H can be obtained if the signif- 
icant wave height is small compared to the water depth and the wave steepness 

is low. Individual wave crests with elevations greater than 1.5 times the 
significant wave height were observed in a number of the laboratory experiments, 
where H,/d <0.3 and d/(gT5) <0.005 for m= 0.0. 


Similar values of maximum wave crest elevation were found for irregular 
waves shoaling on a 1 on 30 laboratory smooth slope. An examination of a typi- 
cal set of results showed that crest elevation increases as the wave shoals, 
then decreases as the wave breaks (Fig. 17). Note that in this figure both the 
crest elevation and the stillwater depth have been normalized by the deepwater 
equivalent significant wave height, Hg. Typical values of the single highest 
crest observed during an irregular wave test of 260 seconds for model time are 
presented in Table 7. Water level time histories at several gages at various 
depths along the 1 on 30 slope are given in Figure 18 for an irregular wave 
condition. 


Jahns and Wheeler (1973) suggested a method for predicting irregular wave 
crest statistics, using the dimensionless parameter (n,-)/Hg. A Rayleigh-type 
distribution is assumed, then an empirical correction factor is recommended 
to increase the probability of occurrence of crests. One of the disadvantages 
of this approach is that the data often dramatically differ from a Rayleigh 
distribution because actual probabilities are orders of magnitude larger than 
predicted. For example, Figure 19 compares crest height distributions for the 
wave records A, B, and C shown in Figure 16 with the Rayleigh distribution. 


The approach taken in this report is to predict crest elevations in terms 
of local parameters for a given probability level in terms of H,, d, and 
T,. Figure 20 shows values of Ne/d for wave records A, B, and C illustrated 
in Figure 16. 


The suggested method for predicting irregular wave crest elevations is to 
use the stream-function results to determine the general form of the equation 
and laboratory results to calibrate the prediction technique for various prob- 
ability levels. The equation used is 


*k 
(ne) p chs Ap eases aN Me an (8) 
H 
Se eons teanh (-a. ta(a; 2)) 


where the parameters Aj and A3 are taken from stream-function results 

(Fig. 8 and Table 2), and A} is an empirical parameter determined from the 
experimental data. The laboratory data (App. C) indicate that A? is a 
function of the probability of exceedance, with A increasing as the prob- 
ability, p, decreases (Fig. 21, Table 8). Equation (8) shows good correla- 
tion with the laboratory data as indicated by the small standard deviation of 
Af about the mean for values of p = 0.01 and greater (Tables 8 and 9). A} at 
p = 0.005 has poor resolution because the data collection runs had only a few 
crests more than 200 crests per run. Figure 22 gives observed and predicted 
crest elevations at the 2-percent probability level of exceedance. 
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Figure 17. Variation in wave crest elevation for a 1 on 30 slope. 


Table 7. Selected irregular wave crest data for waves 
breaking on a 1 on 30 slope. 


Hg/Lo Wave and crest data at the occurrence of the 
highest single crest 


Gone | Mie |G 


0.040 
0.032 
0.024 
0.021 
0.0115 
0.0050 
0.0033 
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Figure 18. Sample of water level time histories for irregular waves at 
selected locations along a 1 on 30 slope. 
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Figure 19. Sample of irregular wave crest elevation distributions 
(from wave records shown in Fig. 16). 
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Figure 20. Sample of wave crest distributions (from wave records 
shown in Fig. 16). 
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Figure 21. Values of Ax for probability of exceedance for 
irregular waves. 
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Table 8. Values of A; for levels of probability 
of exceedance for irregular waves 


(m = 0.0), derived from 122 laboratory 
experiments. 


Probability Coefficient 
of of 


exceedance, p variation 


Table 9. Values of Al for levels of probability 
of exceedance for irregular waves on a 
1 on 30 laboratory slope (m = 0.033), 
derived from 192 laboratory measurements. 
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Figure 22. Observed and predicted crest elevations exceeded by 2 percent 
of the wave crests. 
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The wave crest elevation of irregular waves on a 1 on 30 slope may also be 
described in terms of deepwater wave characteristics. An examination of the 
data indicates that deepwater wave steepness and local stillwater depth control 
the magnitude of crest elevation at a given probability level. For example, 
Figure 23 presents the crest elevation exceeded by 13.5 percent of the wave 
crests for the deepwater wave steepness range 0.0035 < Hj/L, < 0.0075 for 
m = 0.0333 (a 1 on 30 slope with the experimental setup shown in Fig. 6). The 
irregular wave data show a consistent trend with crest elevation increasing in 
the shoaling region and decreasing after breaking. 


d/Ho 


Figure 23. Wave crest elevation exceeded by 13.5 percent of the wave crests. 


2. Crest Elevation Prediction Aids. 


Irregular wave crest elevations may be predicted in terms of local wave 
conditions using equation (8). Figure 24 gives crest elevation exceeded by 
2 percent of the wave crests (p = 0.02) normalized by the water depth. Curves 
of constant H,/d are plotted with d/(gT2) on the abscissa. As a first 
estimate the upper limit of wave height is assumed to be given by equation (5) 
and the crest elevation of this condition is shown as a dashline. Figure 25 
presents crest elevations exceeded by 13.5 percent of the wave crests. Figures 
24 and 25 may be used for waves traveling over a flat bottom and for plane 
slopes as steep as 1 on 30. 


Crest elevations are presented in terms of deepwater wave characteristics 
for incident waves normal to a 1 on 30 slope for the 2- and 13.5-percent 
probability of exceedance levels in Figure 26. These curves were determined 
using the laboratory data presented in Appendix D. A dashline is used to 
indicate a region of uncertainty. 
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Figure 24. Crest elevation at the 2-percent probability level of 
exceedance for irregular waves. 
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Figure 25. Crest elevation at the 13.5-percent probability level of 
exceedance for irregular waves. 
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Figure 26. Wave crest elevations along a 1 on 30 slope for 2- and 13.5-percent 
probability of exceedance levels. 
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3. The Duration of Irregular Wave Crests. 
E | 


The larger waves in an irregular wave train have the greatest nonlinearity, 
so the duration of the crests of these waves tends to be small relative to wave 
period. Figure 27 illustrates this trend by presenting the mean relative crest 
duration for various probability levels of exceedance for waves traveling over 
a flat bottom. For example, 5 percent of the largest crests have a mean dimens tom 
less crest duration T,/T = 0.406, while the average value for all of the wave 
crests (at the 100-percent level) have aba//ie = 0.473. 


VI. COMPARISON OF RESULTS FOR PROTOTYPE DATA 


The techniques presented in the previous section are to be used for the 
future prediction of wave characteristics in field situation. In this section, 
field data collected at the FRF in Duck, North Carolina, are analyzed in the 
same manner as the laboratory data previously described and comparisons are 
presented to demonstrate the applicability of the predictive techniques. 


1. Field Conditions. 


Four Waverider buoys were used to obtain the prototype data in an area 
where the hottom contours are essentially straight and parallel. Table 10 gives 
the locations of the gages and the mean water depths. The data were collected 
in the third week of October 1980 during a severe storm that produced wave 
heights exceeding 4 meters. The data set has widely varying significant wave 
height and period of peak energy density combinations resulting as the storm 
moved into and through the study area. Table 10 summarizes the range of wave 
conditions used in this report. 


2. Prototype Verification. 


The field data were analyzed in the same manner as the laboratory data used 
in formulating the predictive parameters. The observed crest elevations were 
compared with those predicted by equation (8), using values of AT for m= 0.0 
(Table 8), since the bottom slope was very gentle. Table 11 shows the results 
of this comparison for decreasing exceedance probabilities. It should be noted 
that the predicted crest elevations are slightly conservative, with the coef- 
ficient of variation approximately 8.5 percent. Figures 28 and 29 illustrate 
the predictive ability of equation (8); the diagonal line indicates perfect 
agreement. The scatter decreases for higher probability of exceedance levels, 
in part, because the data for higher probability levels are more statistically 
reliable. These figures also show that the scatter from the prediction line 
increases slightly as H, /d increases. This increase in conservatism with 
increasing wave height-to-depth ratio may be due to the fact that the waves 
are interacting with a bottom that is slightly sloping and comprised of a 
porous sandy material (rather than a rigid floor used in the laboratory 
experiments). The directional characteristics of the prototype waves could 
also be influencing these results because a spectrum of wave directions was 
present, even though the primary wave direction was along the line of the 
Waveriders. 
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Percentage of highest wave crests 
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Dimensionless wave crest duration for irregular waves (m = 0.0). 


Figure 27. 


Table 10. Summary of prototype wave conditions. 


Waverider 

distance 

from shore 
(km) 


0.0016 
0.0016 
0.0013 
0.0013 


Table 11. Values of crest heights observed and predicted 
for waves at Duck, North Carolina. ! 


Probability - Mean ratio of Ai Coefficient 
of observed and Observed of 
exceedance predicted variation 

crest elevations (pet) 


lyalues derived from 100 20-minute wave records taken in 
October 1980. 


* 
2values of A, from Table 8 used to make predictions. 
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VII. SUMMARY AND CONCLUSIONS 


The stream-function tables of Dean (1974) are used to develop empirical 
formulas for estimating monochromatic wave crest elevations and durations. 
The resulting formulas are calibrated with laboratory data. The method for 
predicting crest elevations was found to work well for waves traveling over 
a flat bottom or at the breaker point on plane smooth slopes. The equation 
for predicting crest elevations was calibrated for irregular wave conditions 
traveling over a flat bottom and shoaling on a 1 on 30 plane slope. The 
irregular wave data taken with the 1 on 30 slope were also expressed in terms 
of deepwater conditions and it was found that deepwater steepness and local 
water depth could be used to predict crest elevations. 


Prediction equations developed from two-dimensional irregular laboratory 
waves were found to be conservative in predicting crest elevations when 
compared to field measurements taken at CERC's Field Research Facility, Duck, 
North Carolina. The conservatism is slight and increases with increasing 
H,/d. 
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APPENDIX A 


MONOCHROMATIC TESTS RESULTS 


BREAKING MONOCHROMATIC WAVE DATA 
COT(M)= 95 
DATA BY SINGAMSETTI AND WIND(1980) 


HB(CM) T(SEC) DB(CM) DB/LO HB/LO ETA C/H TC/T 
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9.6 1.2 10.2 0.040 0. 0376 0. 608 0.453 
a8) io Be) 7.2 0.019 0. 0252 0.654 90.412 
10.5 1.04 3.0 0. 048 0. 0628 0.653 0.373 
ee G? 1. 28 13.0 0.051 0. 0621 QO. 645 0. 407 
9.0 1.04 7.7 0. 046 0. 0538 0. 639 0. 408 
Ais 7 1. 04 9.6 0.057 0. 0696 0.675 0.419 
18.3 1.72 Was 0.042 0. 0398 0.581 0.392 
2.4 We Ve DES 0.024 0. 0270 0. 628 0. 42 
14.2 1.2 13.6 0.053 0. 0631 0. 648 0. 397 
2.1 1.28 10.1 0. 039 0. 0471 0. 674 0. 440 
8.9 1.28 8.0 0.031 0. 0347 0. 657 0. 445 
9.2 1.55 9.6 0.025 0. 0246 0.619 0.402 
8.7 1.04 8.0 0. 047 0.0518 0. 638 0. 398 
15.0 1. 28 13.0 0.051 0. 0582 0.615 0. 370 
Gat 1.04 Pe 0 0. 058 9. 0453 0. 635 0.43 
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No EI he WE Mig @) 0. 038 0. 0268 0. 586 0.498 

SUM, SUI11, NO= 9.4720 1 0858 @1 

SD,SD1 =9. 036970. 09162 

TC/T N, AVERAGE SD= 21 1.8236 0. 3165 

BREAKING MONOCHROMATIC WAVE DATA 
cOT(M)= 10 
DATA BY SINGAMSETTI AND WIND( 1980) 
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SUM, GUIML, WO= 0.8181 1.2420 22 
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APPENDIX B 


IRREGULAR WAVE CREST ELEVATIONS 
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APPENDIX E 


CNOIDAL WAVE TESTS RESULTS 


Wave Wave Water Relative 
height period depth crest 
height 
(cm) (s) (cm) n/H 
11.278 2.30 44.05 0.6260 
F583 2.80 44.05 0.6658 
11.433 3.40 44.05 0.7526 
11.468 4.10 44.05 0.7678 
11.590 4.80 44.05 0.7767 
11.464 5.70 44.05 0.8118 
Za © 1.60 44.05 0.6330 
Pale @7/ 2.00 44.05 0.6470 
23.38 2.40 44.05 0.6895 
ZBI. 2.40 44.05 0.7258 
21.88 SQ) 44.05 0.7724 
23.104 Soo) 44.05 0.7228 
ihe as) 4.00 44.05 0.8036 
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APPENDIX F 


ESTIMATING THE WAVE CREST ELEVATION AT BREAKING ON A PLANE SLOPE 


The equations of Singamsetti and Wind (1980) can be used with the methods 
developed in this report to predict the crest elevation at the breaking point 
for monochromatic waves normally incident to a slope. This technique uses as 
input water depth, d, wave height, H, and period, T, in a constant depth 
region offshore of a plane slope, m. Deepwater wave height is found from 


H 
He = ——————— (F-1) 


1.0 


eR MERIC HLL, Neate elen 
21d 2 
are 0 i Ga 


where Ly is the wavelength from Airy theory 
2 
2 ait 21d Ae 
iLyf\ On tanh i (F-2) 


Breaking wave height may be determined from 


0.031 -0.254 
Hy = HS 0.575(m) (Hi/L,) (F-3) 
where L, is deepwater Airy wavelength 
2 
ge = 
Lois an (F-4) 
Depth at breaking, dp, is found from 
-0.155 0,130 
d, = H, (1.067) (@m) CHg/ES) (F-5) 
fo a m a 
2 BG 5 


-0.107 0.237 
= H, (1.761) (m) (aS /ie) 


a 
ion 
I 


for m =< 


The value of dp/L, is used to determine A, and A from Figure 6 in the 
test of this report; A, is determined from the value of the beach slope 
using Figure 7. Crest height at breaking, (n.)p, is then given by 


A, H, 


(F-6) 
1.0 + tanh [ote Ln (1s | 
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